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Abstract 

We calculate the nuclear energy density functional relevant for N=Z even-even nuclei in 
the systematic framework of chiral perturbation theory. The calculation includes the one- 
pion exchange Fock diagram and the iterated one-pion exchange Hartree and Fock diagrams. 
From these few leading order contributions in the small momentum expansion one obtains al- 
ready a very good equation of state of isospin symmetric nuclear matter. We find that in the 
region below nuclear matter saturation density the effective nucleon mass M*(p) deviates by 
at most 15% from its free space value M, with 0.89M < M*{p) < M for p < 0.11 fm~ 3 and 
M*(p) > M for higher densities. The parameterfree strength of the (V/9) 2 -term, F\?(kf), is 
at saturation density comparable to that of phenomenological Skyrme forces. The magnitude 
of Fjikf) accompanying the squared spin-orbit density J 2 comes out somewhat larger. The 
strength of the nuclear spin-orbit interaction, F so (kj), as given by iterated one-pion exchange 
is about half as large as the corresponding empirical value, however, with the wrong negative 
sign. The novel density dependencies of M*(p) and Fv iSO ,j(fc/) as predicted by our param- 
eterfree calculation should be examined in nuclear structure calculations (after introducing 
an additional short range spin-orbit contribution constant in density). 

PACS: 12.38.Bx, 21.30.Fe, 21.60.-n, 31.15.Ew 
Keywords: Nuclear energy density functional; Density-matrix expansion; Perturbative chiral pion- 
nucleon dynamics 



1 Introduction 

Among the various phenomenological interactions that have been used extensively in the de- 
scription of nuclei, the Skyrme force H has gained much popularity because of its analytical 
simplicity and its ability to reproduce nuclear properties over the whole periodic table within the 
self-consistent Hartree- Fock approximation. Several Skyrme parameterizations have been tailored 
to account for single particle spectra [[J , giant monopole resonances || or fission barriers of heavy 
nuclei ||. Recently, a new Skyrme force which also reproduces the equation of state of pure 
neutron matter up to neutron star densities, p n ~ 1.5 fm -3 , has been proposed in ref.[|5| for the 
study of nuclei far from stability. A microscopic interpretation of the various parameters entering 
the effective Skyrme forces is generally put aside. Sometimes the energy density functional is 
just parameterized without reference to any effective (zero-range) two-body interaction in order 
to avoid the complete fixing of time-reversal-odd terms by the Pauli-exclusion principle |J. 

!Work supported in part by BMBF, GSI and DFG. 
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Another widely and successfully used approach to nuclear structure calculations are relativistic 
mean-field models [0. In these models the nucleus is described as a collection of independent 
Dirac-particles moving in self-consistently generated scalar and vector mean-fields. The footprints 
of relativity become visible through the large nuclear spin-orbit interaction which emerges in that 
framework naturally from the interplay of the two strong and counteracting (scalar and vector) 
mean-fields. The corresponding many-body calculations are usually carried out in the Hartree 
approximation, ignoring the negative-energy Dirac-sea. For a recent review on self-consistent 
mean-field models for nuclear structure, see ref.Q. In that article the relationship between the 
relativistic mean-field models and the Skyrme phenomenology is also discussed. 

The first conditions to be fulfilled by any phenomenological nucleon-nucleon interaction come 
from the (few empirically known) properties of infinite nuclear matter. These are the saturation 
density p = 2/c 3 /37r 2 , the binding energy per particle —E(kf ) and the compression modulus 
K = kj Q E"(kfo) of isospin symmetric nuclear matter as well as the asymmetry energy A(kfo). 
In general, Skyrme forces involve several more parameters related to terms in the energy density 
functional which vanish identically in homogeneous nuclear matter (like the spin-orbit coupling 
proportional to the density- gradient). 

In a recent work ||, we have used chiral perturbation theory for a systematic treatment of 
the nuclear matter many-body problem. In this calculation the contributions to the energy per 
particle, E(kf), originate exclusively from one- and two-pion exchange between nucleons and they 
are ordered in powers of the Fermi momentum kf (modulo functions of kf/m n ). It has been 
demonstrated in ref. || that the empirical saturation point (po — 0.17 fm -3 , E(kf ) ~ —15.3 MeV) 
and the nuclear matter compressibility K ~ 255 MeV can be well reproduced at order 0(kf) in 
the small momentum expansion with just one single momentum cut-off scale of A ~ 0.65 GeV 
which parameterizes the necessary short range NN-dynamics. Most surprisingly, the prediction 
for the asymmetry energy, A(kfo) = 33.8 MeV, is in very good agreement with its empirical 
value. Furthermore, as a nontrivial fact, pure neutron matter is predicted to be unbound and 
the corresponding equation of state E n {k n ) agrees roughly with that of sophisticated many-body 
calculations for low neutron densities p n < 0.25 fm -3 . In a subsequent work ||, the momentum 
and density dependent (real) single-particle potential U(p,kf) has been calculated in the same 
framework. It was found that chiral In- and 27r-exchange give rise to a potential depth for a 
nucleon at the bottom of the Fermi sea of {7(0, fc/o) = —53.2 MeV. This value is in very good 
agreement with the depth of the empirical optical model potential and the nuclear shell model 
potential. In fact very similar results for nuclear matter can be obtained already at order O(kj) 
in the small momentum expansion (by dropping the relativistic l/M 2 -correction to l7r-exchange 
and the irreducible 27r-exchange of order 0(k^)) with a somewhat reduced cut-off scale of A ~ 
0.61 GeV. Detailed results of this 0(/c|)-calculation can be found in ref. |TIJ (it is named "mean- 
field treatment of the NN-contact interaction" therein). 

Given the fact that many properties of nuclear matter can be well described by chiral itN- 
dynamics treated perturbatively up to three-loop order it is natural to consider in a further step 
the energy density functional relevant for inhomogeneous many- nucleon systems (i.e. finite nuclei). 
We will restrict ourselves here to the isospin symmetric case of equal proton and neutron number 
N=Z. The aim of this work is to calculate density-dependent generalizations of those combinations 
of Skyrme parameters which belong to terms specific for the inhomogeneous many-nucleon system, 
such as (Vp) 2 . This novel density-dependence is a consequence of the finite range character of the 
17T- and 27r-exchange interaction. We stress already here that our results for the density- dependent 
strength functions are completely parameterfree. In particular, they are independent of the cut-off 
scale A adjusted in ref. [EL [TTJ to the binding energy per particle —E(kf ) = 15.3 MeV. 



Our paper is organized as follows. In the next section we recall the essential results of the 



density-matrix expansion of Negele and Vautherin [TT| which provides the adequate technical 
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framework to compute the nuclear energy density functional. In section 3, we then present ana- 
lytical formulas for the various strength functions F T (kf), F d (kf), F so (kf) and Fj(kf) entering the 
nuclear energy density functional. These expressions are derived (exclusively) from the two-loop 
l7r-exchange Fock diagram and the three-loop iterated l7r-exchange Hartree and Fock diagrams. 
Section 4 is devoted to a discussion of our results and finally section 5 ends with a summary and 
an outlook. 



2 Density-matrix expansion and energy density functional 

The starting point for the construction of an explicit nuclear energy density functional is the 
density-matrix as given by a sum over the occupied energy eigenfunctions \l/ a of this many-fermion 



system. According to Negele and Vautherin [1 1J the bilocal density-matrix can be expanded in 
relative and center-of-mass coordinates, a and r, as follows: 
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where the functions ji(akf) are ordinary spherical Bessel functions. The other quantities appearing 
on the right hand side of eq.(l) are the local nucleon density: 
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written here in terms of the local Fermi- momentum kf(f), the local kinetic energy density: 

E Wl(f)-VvMr), (3) 
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and the local spin-orbit density: 



J(f) = E Vl(r)ie x W a (f) 



(4) 
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For notational simplicity we have dropped their argument r in eq.(l) and will do so in the following. 
It is important to note that a pairwise filling of time-reversed orbitals a has been assumed in eq.(l). 
If the many-body ground state is not time-reversal invariant (as it is the case for odd nuclei and 
for rotating nuclei) various additional time- reversal-odd fields come into play ||. The local spin- 
orbit density J(r ) is non-zero for spin-unsaturated shells. In such a situation the density-matrix 
is no longer a scalar in spin-space but has also a vector part. The Fourier transform of the 
density-matrix eq.(l) with respect to both coordinates a and f defines the medium insertion for 
the inhomogeneous many- nucleon system characterized by the time-reversal-even fields p(r), r(f) 
and J(f ): 
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The double line in the left picture of Fig. 1 symbolizes this medium insertion together with the 
assignment of the out- and in-going nucleon momenta p ± q/2. The momentum transfer q is 



provided by the Fourier components of the inhomogeneous (matter) distributions: p(r), r(r ) and 
J(r). As a check one verifies that the Fourier transform (1/2-7T 3 ) J d 3 pe~ tp ' a of the (partly very 
singular) expression in the curly brackets in eq. (5) gives exactly the right hand side of the density- 
matrix expansion in eq.(l). For homogeneous nuclear matter (where r = 3pk 2 /5 and Vp = J = 0) 
only the familiar step-function 8(kf — |p |) remains from the medium insertion eq.(5) as the density 
of nucleon states in momentum space. 

Going up to second order in spatial gradients (i.e. deviations from homogeneity) the energy 



density functional relevant for N=Z even-even nuclei reads |T2 



S[p,r,J] = pE(k f ) + 



\P k ) 



5k 2 , s 

' +Fr(kf) 



2M 56M 3 

+(Vp) 2 F v (k f ) + Vp ■ JF so (k f ) + J 2 Fj(k f ) . (6) 

Here, E[kf) is the energy per particle of isospin symmetric nuclear matter evaluated at the local 
Fermi momentum kf(f). The (small) relativistic correction term — 5k 2 /56M 3 has been included 
in eq.(6) for the following reason. When multiplied with —3pk 2 /5 it cancels together with the fore- 
going term 1/2M the relativistically improved kinetic energy in E(kf) (see eq.(5) in ref.pj). The 
functions F T (kf), F^(kf), F so (kf) and Fj(kf) arising from NN-interactions encode new dynami- 
cal information specific for inhomogeneous many-nucleon systems. In Skyrme parameterizations 
F T {kf) depends linearly on the (local) density p = 2k 3 /3n 2 whereas Fv iSO ,j(^/) are just constants. 
Note that F so {kf) gives the strength of the nuclear spin-orbit coupling while F\r(kf) is responsible 
for the formation of the nuclear surface. Variation of the energy density functional S[p, r, J] with 
respect to single-particle wavefunctions under the condition that these are normalized to unity 



leads to self-consistent density dependent Hartree-Fock equations [|12 |. 

Returning to the medium insertion in eq.(5) one sees that the strength function F T {kf) emerges 
via a perturbation on top of the density of states 0{kf — \p\). The single particle potential in 
nuclear matter can actually be constructed in the same way by introducing a delta-function like 
perturbation [[|. Consequently, the strength function F T (kf) can be directly expressed in terms 
of the real part U(p, kf) of the momentum and density dependent single particle potential as: 

FAkf) = w f l d pp 2 ( 5 p 2 - 3k2 f) u (p> k f)- ( 7 ) 

In eq.(5) the term r — 3pk 2 /5 is accompanied by — V 2 p/4. Performing a partial integration of 
the energy / d 3 rS one sees immediately that part of the strength function F\/(kf) is given by the 
p-derivative of F r (A;/)/4. These considerations lead to the following decomposition: 

where Fd(kf) comprises all those contributions for which the (Vp) 2 -factor originates directly from 
the interactions. An example for this mechanism will be explained in the next section. 

As a check on the present formalism (summarized in eq.(5)) we rederived the Skyrme energy 
density functional (eq.(5.87) in ref. |12|) from the matrix elements of the underlying two-body 
potential (eq.(4.105) in ref.|l2|) in a purely diagrammatic framework. In the next section we 



use the same formalism to compute the nuclear energy density functional eq.(6) from one- and 
two-pion exchange diagrams. 
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3 Diagrammatic calculation 



In this section we present analytical formulas for the four strength functions F T (kf), Fa(kf), 
F so (kf) and Fj(kf) as derived from the two-loop one-pion exchange Fock diagram and the three- 
loop iterated one-pion exchange Hartree and Fock diagrams. These graphs are shown in Fig. 1. 
We give for each diagram only the final result omitting all technical details related to extensive 
algebraic manipulations and solving elementary integrals. 




( > f + a/2 



Fig. 1: Left: The double line symbolizes the medium insertion defined by eq.(5). Next are shown: 
The two-loop one-pion exchange Fock-diagram and the three-loop iterated one-pion exchange Hartree- 
and Fock- diagrams. The combinatoric factors of these diagrams are 1/2, 1/4 and 1/4, in the order 
shown. Their isospin factors for isospin symmetric nuclear matter are 6, 12 and —6, respectively. 



3.1 One-pion exchange Fock diagram with two medium insertions 

The non-vanishing contributions read: 

Fr(kf) = 70^5 + 2W " 1 " 20 W arctan2 M + (J - 6m 2 + ±) ln(l + An")) , (9) 



g\ (10 + 24m 2 1 Wi , 2 



(16tt/ w 

FAkf) = (TOi f^? + ^ ln(1 + 4M2) i ' (10) 
with u = kf/m n the ratio of the two small scales inherent to our calculation. The expression for 
F T (kf) in eq.(9) follows simply from inserting the static l7r-exchange single particle potential 
(eq.(8) in ref. |]] in the limit M — * oo) into the "master formula" eq.(7). The vanishing of 
Fd(kf) has the following reason. The momentum transfer ±q at the upper and lower medium 
insertion does not flow into the exchanged virtual pion line (because of momentum conservation 
at each interaction vertex). Therefore there is no factor of q 2 which could produce via Fourier 
transformations a (Vp) 2 -factor. The spin-orbit strength F so (kf) vanishes as a result of the spin- 
trace: tr[<r • (pi - p 2 ) a ■ (p lj2 x J ) a ■ (pi - p 2 )] = 0. 



3.2 Iterated one-pion exchange Hartree diagram with two medium 
insertions 

We find the following closed form expressions: 

Fr(kf) = /^iJ 1 ^ - (350 + W) arctan 2u 
(Stt) 3 ^/^) 4 I 3 

+444m- — + — - + - -« ln l + 4« 2 , 11 

4m \16m 3 8m 2 / v ') v ' 



= ^m\u " cten2u - is? ln(1 + 4u > - i2 ( i + tf)4 ■ < 12 > 

Again, F T (kf) in eq.(ll) stems from inserting the two-body potential ^(p, fc/) (eq.(9) in ref.|J) 
into the "master formula" eq.(7). Note that any p-independent contribution, in particular the 
cut-off dependent term eq.(17) in ref.||, drops out. The vanishing of Fj(kf) results from the 
spin-trace over a nucleon ring being equal to zero (as demonstrated at the end of section 3.1). Let 
us briefly explain the mechanism which generates the strength function Fd(kf). The exchanged 
pion-pair transfers the momentum q between the left and the right nucleon ring. This momentum 
q enters both the pseudovector 7riV-interaction vertices and the pion propagators. After expanding 
the inner loop integral to order q 2 the Fourier transformation in eq.(5) converts this factor q 2 
into a factor (Vkf) 2 . The rest is a solvable integral over the product of two Fermi surfaces. The 
spin-orbit strength F so (kf) arises from the spin-trace: tx[a ■ (I + q/2) a ■ (I — q/2) a • (pi )2 x J)] = 
2i(qxl)- (pip x J ) where q gets again converted to Vkf by Fourier transformation. The remainder 
is a solvable integral over delta-functions and derivatives thereof. 



3.3 Iterated one-pion exchange Fock diagram with two medium inser- 
tions 

We find the following contributions from the last diagram in Fig. 1 with two medium insertions at 
non-neighboring nucleon propagators: 
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The basic mechanisms which lead to these results are the same as explained before. Concerning 
kinematics and spin-algebra the iterated l7r-exchange Fock diagram is somewhat more involved 
than the Hartree diagram. Even though all occurring inner c? 3 /-loop integrals can be solved in 
closed form there remain some non-elementary integrals from the integration over the product of 
two Fermi spheres of radius kf. 



3.4 Iterated one-pion exchange Hartree diagram with three medium 
insertions 



In this case we find the following contributions: 
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with the auxiliary functions H = \n(u + xy) — \n(u — xy) and s = xy + v^^ 2 — + x 2 y 2 . The 
quantity s has the geometrical meaning of the distance between a point on a sphere of radius u and 
an interior point displaced at a distance x from the center of the sphere. In the same geometrical 
picture y denotes a directional cosine. 
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Here we have introduced the partial derivatives s' = uds/du and s" = u d s/du 
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The last contribution Fj(kf) in eq.(21) is obtained when both insertions proportional to cr-(pi i2 x J ) 
(producing, after integration, the overall J 2 -factor) are under a single spin-trace. For the other 
two possible combinations the spin-traces are equal to zero. 



3.5 Iterated one-pion exchange Fock diagram with three medium in- 
sertions 

The evaluation of this diagram is most tedious. It is advisable to split the contributions to the four 
strength functions F T (kf), F d (kf), F so {kf) and Fj(kf) into "factorizable" and " non-factorizable" 
parts. These two pieces are distinguished by whether the nucleon propagator in the denominator 
can be canceled or not by terms from the product of 7riV-interaction vertices in the numerator. 
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We find the following "factorizable" contributions 
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The " non-factorizable" contributions read on the other hand: 
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(30) 



with the auxiliary function t = xz + Vw 2 — x 2 + x 2 z 2 and its partial derivatives t' = udt/du and 
t" = u 2 d 2 t/du 2 . For the numerical evaluation of the ofa/iiz-double integrals in eqs. (27-30) it is 
advantageous to first antisymmetrize the integrands both in y and z and then to substitute z = 
Vy 2 C 2 + 1 — V 2 - This way the integration region becomes equal to the unit-square < y, ( < 1. 

Obviously, the analytical results presented in this section do not involve any adjustable param- 
eter. Only well-known physical quantities like the nucleon axial- vector coupling constant qa = 1-3, 
the nucleon mass M = 939 MeV, the pion decay constant /„- = 92.4 MeV and the (neutral) pion 
mass m n = 135 MeV enter. 

Let us end this section with general power counting considerations for the nuclear energy 

— * 

density functional £[p,r,J]. Counting the Fermi momentum kf, the pion mass m T and a spatial 
gradient V collectively as small momenta one deduces from eqs. (2,3,4) that the nucleon density 

— # 

p(r), the kinetic energy density r(f) and the spin-orbit density J(r) are quantities of third, fifth 
and fourth order in small momenta, respectively. With these counting rules the contributions 
from l7r-exchange to the nuclear energy density functional S[p,r,J] are of sixth order in small 
momenta while all contributions from iterated l7r-exchange are of seventh order. Concerning Non- 
interactions induced by pion-exchange the nuclear energy density functional presented here is in 
fact complete up-to-and-including seventh order in small momenta. 



4 Results and discussion 



In this section we present and discuss our numerical results using the input parameters just 
mentioned. Returning to the energy density functional S[p,r, J] in eq.(6) one observes that the 
expression in square brackets multiplying the kinetic energy density r(r ) has the interpretation 
of a reciprocal density dependent effective nucleon mass: 



M*(p) = M 



5*? 
28M 2 



+ 2MF T (k f ) 



i -l 



(31) 
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We note as an aside that this effective nucleon mass M*(p) is conceptually different from the so- 
called "Landau" mass which derives from the slope of the single particle potential U(p, kf) at the 
Fermi surface p = kf. Only if the (real) single particle potential has a simple quadratic dependence 
on the nucleon momentum, U(p, kf) = Uo(kf) + Ui(kf)p 2 , these two variants of effective nucleon 
mass agree with each other. 
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Fig. 2: The effective nucleon mass M*(p) divided by the free nucleon mass M versus the nucleon 
density p. The dotted line corresponds to the fit: M*(p)/M = 1 — 3.054 fm 2 ■ p 2//3 + 6.345 fm 3 ■ p. 

In Fig. 2 we show the ratio effective over free nucleon mass M*(p)/M as a function of the 
nucleon density p = 2A; 3 /37r 2 . One observes a reduced effective nucleon mass 0.89M < M*(p) < M 

for densities p < 0.11 fm -3 and an enhanced effective nucleon mass M*(p) > M for higher densities. 
In the region below the nuclear matter saturation density p < po = 0.174fm~ 3 relevant for nuclear 
structure the deviations of the effective nucleon mass M*(p) from its free space value M do no 
exceed ±15%. Let us give a qualitative explanation for the (somewhat unusual) behavior of 
the curve in Fig. 2. Consider the iterated l7r-exchange Hartree diagram in Fig. 1 at sufficiently 
high densities such that the pion mass m v can be neglected against the Fermi momentum kf. 
In this (limiting) case the 7riV-interaction vertices get cancelled by the pion propagators. One 
is effectively dealing with a zero-range contact interaction in second order which according to 
Galitskii's calculation from 1958 [[13], Q generates an enhanced in-medium mass. In this sense 
the curve in Fig. 2 delineates the two density regimes kf < \[?>m^ and kf > ^/3m^ where the 
(qualitative) behavior in the latter is ruled by Galitskii's second order result. Interestingly, a 



recent large scale fit of 1888 nuclide masses by Pearson et al. ]15|] using a " Hartree- Fock nuclear 
mass formula" has given an effective nucleon mass of M*(po) = 1.05M. This value is comparable 
with our parameterfree result M*(po) = 1.15M. 
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Fig'. 3: The strength function F^(kf) related to the (V p) 2 -term in the nuclear energy density 
functional versus the nucleon density p = 2kj/3ir 2 . An accurate fit of the full line is: Fy(kf) = 
45.43 Me Vfm 4 • p -1 / 3 — 0.229 Me Vfm 2 ■ p -1 . The three horizontal dashed lines show the constant 
values F v (k f ) = [9h - (5 + 4x 2 )t 2 ]/64 of the Skyrme forces Sly §J, SIII /@/ and MSk gBjj. 

The dotted line in Fig. 2 corresponds to the fit: M*(p)/M = 1 - 3.054 fm 2 • p 2 / 3 + 6.345 fm 3 ■ 
p, which may be useful for applications in nuclear structure calculations. In this context we 
mention also the fitted form of the underlying nuclear matter equation of state fllUfl : E(kf) = 
111.63 Me Vfm 2 ■ p 2 / 3 - 752.82 Me Vfm 3 • p + 832.74 Me Vfm 4 ■ p 4 / 3 . 

Next, we show in Fig. 3 by the full line the strength function Fy(kf) belonging to the (Vp) 2 - 
term in the nuclear energy density functional eq.(6) versus the nucleon density p = 2/c 3 /37r 2 . The 
three horizontal dashed lines represent the constant values F v (/c/) = [9ti — (5 + 4x 2 )t 2 ]/64 of the 
Skyrme forces Sly ||, SIII and MSk fll5 |. In the case of Sly and MSk we have performed averages 
over the various parameter sets Sly4-7 and MSkl-6. At nuclear matter saturation density po = 
0.174 fm -3 our parameterfree prediction Fv(/c/o) = 80.1 Me Vfm 5 is comparable to these empirical 
values. The strong increase of the strength function Fy(/c/) with decreasing density has to do 
with the presence of a small mass scale, m n = 135 MeV, and with associated chiral singularities 
(of the form m~ 2 and m' 1 ). We will come back to this issue again towards the end of this section. 
An accurate fit of the full line in Fig. 3 is: F v (k f ) = 45.43 Me Vfm 4 • p^ 1 / 3 - 0.229 Me Vfm 2 ■ p~ l . 
It may be useful for applications in nuclear structure calculations. Note also that the relevant 
contribution to the central sing le-particle potential, -(Vp) 2 dF v (k f )/dp - 2V 2 pF v (k f ), receives 
only little weight from very low densities. Therefore the deviation of the strength function Fy(kf) 
from a constant may be less dramatic in practice than it appears on first sight from Fig. 3. 
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Fig. 4: The strength function F so (kf) related to the spin-orbit coupling term in the nuclear energy 
density functional versus the nucleon density p = 2A; 3 /37r 2 . The dotted line corresponds to the fit: 
F so (k f ) = 1.898 Me Vfm 3 • p~ 2 / 3 - 29.37 Me Vfm 4 • p" 1 / 3 . The three horizontal dashed lines show 
the constant values F s0 (kf) = 3Wo/4 of the Skyrme forces Sly jj^j, SIII ||/ and MSk |7^/ . The 
dashed-dotted line shows the contribution from irreducible lis-exchange written in eg. (33) for a 
cut-off A = 0.65 GeV. 



The full line in Fig. 4 shows the result of iterated l7r-exchange for the strength function F so (kf) 
belonging to the spin-orbit coupling term in the nuclear energy density functional. For comparison 
we have drawn the constant values F so (kf) = 3^0/4 of the three Skyrme forces Sly SIII @ 
and MSk jTS] (horizontal dashed lines). One observes that the strength of the nuclear spin-orbit 
interaction as generated by iterated l7r-exchange is at po about half as large as the corresponding 
empirical value, however, with the wrong negative sign. This "negative" result is dominated by 
the contribution of the iterated l7r-exchange Hartree diagram with two medium insertions (see 
eq.(13)). For example, one obtains numerically from eq.(13) at saturation density po = 0.174fm -3 
(where u = kf /m n = 2.0) the negative value F so (kf ) = —83.7 Me Vfm 5 . The other diagrams 
with lower spin- and isospin weight factors reduce this number approximately half in magnitude. 
The "negative" result for F so (kf) is to some extent already indicated by the calculation of the 
momentum and density dependent nuclear spin-orbit strength Ui s (p,kf) in ref. |16|- Going back 
to the medium insertion in eq.(5) one learns that only the values of Ui s (p,kf) near the Fermi 
surface p = kf will contribute to F so (kf). As a matter of fact the curves in Fig. 7 of ref.fTB" 
drop from positive to negative values when p runs from zero to fc/o = 272.7 MeV. Actually, for 
the contributions to F so (kf) from diagrams with two medium insertions eqs. (13,16) the following 
relationship holds: 



Fso(kf) 



7T 



dU ls (p, kj 
dk f 



+ 



k f d 2 Ui s (p,k f ) 



dpdkj 



p=k f 



(32) 
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to be applied to the expressions {7 Z °' (p, kf) in eqs.(9,17) of ref . fl~6| . 

It is well-known that irreducible two-pion exchange generates (via relativistic 1/M-corrections) 
spin-orbit amplitudes in the T- matrix of elastic nucleon-nucleon scattering [17]. Their effect on the 



nuclear spin-orbit interaction has been calculated in ref. |1| . Inserting the expression (p, k 



in eq.(18) of ref.|18[ into the "master formula" eq.(32) one derives the following contribution from 
irreducible 27r-exchange to the spin-orbit strength function: 
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(33) 



Here, A is a momentum cut-off which has been used to regularize the linear divergences of the 
irreducible 27r-exchange (triangle and box) diagrams. In dimensional regularization (employed 
in eqs. (22,23) of ref.|l7]]) such linear divergences are not visible. The dashed-dotted line in Fig. 4 
shows the relatively small contribution of irreducible 27r-exchange to the spin-orbit strength F so (kf) 
for a cut-off scale of A = 0.65 GeV ||. We note that without the zero-range A-dependent term in 
eq.(33) the dashed-dotted curve would be shifted downward by 45.7MeVfm 5 to negative values. 
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Fig. 5: The strength function Fj(kf) accompanying the squared spin-orbit density J 2 in the nuclear 
energy density functional versus the nucleon density p = 2/c|/37r 2 . An accurate fit of the full line 

1/,s . The three horizontal dashed lines show 



is: Fj(k f ) = 12.80 Me Vfm 7/2 -p 
the constant values Fj(kf) = [ti(l 
and Sly j^j. 



x /2 + 7.041 MeVfm 4 ■ p~ 



2xi) - t 2 (l + 2x 2 )]/32 of the Skyrme forces MSk j^j, SIII 



In this context we mention also the relativistic l/M 2 -correction to F so (kf) from the l7r-exchange 

/) in eq.(6) of ref. [[RJ into the "master formula" 



Fock diagram. Inserting the expression (p, k 
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eq.(32) leads to the simple result: F so (kf) = 5r^[ln(l + 4m 2 ) — 4u 2 ]/(16M/,w) 2 , with u = kf/m n . 
As expected, this contribution with F so (2m 7T ) = — 0.86MeVfm 5 is negligibly small. The full line 
in Fig. 4 is accurately fitted by: F so (k f ) = 1.898 Me Vfm 3 • p~ 2 / 3 - 29.37 Me Vfm 4 • p" 1 / 3 . 

Finally, we show in Fig. 5 the strength function Fj(kf) accompanying the squared spin-orbit 
density J 2 in the nuclear energy density functional versus the nucleon density p = 2/c 3 /37r 2 . For 
comparison we have drawn the constant values Fj(kf) = [ti(l — 2xi) — t 2 (l + 2x 2 )]/32 of the three 
Skyrme forces MSk fl~5fj , SIII and Sly (dashed lines). One observes that our prediction for 



Fj(kf) is considerably larger. Again, there is a strong rise of the strength function Fj(kf) as one 
goes down to very low nucleon densities p < po/10. This time the dominant contribution comes 
from the iterated l7r-exchange Hartree diagram with three medium insertions eq.(21), which gives 
numerically at saturation density Fj(2m,) = 52.5 Me Vfm 5 . It should also be noted that the J 2 - 
term in the energy density functional is often neglected in nuclear structure calculations. The full 
line in Fig. 5 is accurately fitted by: Fj(k f ) = 12.80 Me Vfm 7/2 ■ p- 1 ' 2 + 7.041 Me Vfm 4 -p~ 1/3 , which 
may be useful for implementation into nuclear structure calculations. Note also that the J 2 -term 
in the energy density functional gives rise to an additional spin-orbit single-particle field of the 
form 2Fj(kf)J. According to our calculation this additional spin-orbit field would be rather large 
and strongly density dependent. 

The full curves in Figs. 3,5 show a strong increase as the density p tends to zero. Although not 
visible, each curve approaches a finite value at p = 0. One can analytically derive the following 
low density limits: 
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= 571.3 Me Vfm 5 , (34) 

= 339.2 Me Vfm 5 , (35) 
552.2 Me Vfm 5 , (36) 



^ (0 > = -d^ = - 101 - 4MeVfm5 - (37 > 

to which only the diagrams with two medium insertions contribute. The large numbers in eqs.(34- 
37) arise from negative powers of the pion mass m, (so-called chiral singularities). The most 
singular m~ 2 -terms can be traced back to the l7r-exchange Fock diagram. It is important to keep 
in mind that if pionic degrees of freedom are treated explicitly in the nuclear matter problem the 
low density limit is realized only at extremely low densities kf << m,/2. Often, the opposite limit 
where the pion mass m, can be neglected against the Fermi momentum kf is already applicable at 
the moderate densities relevant for conventional nuclear physics. This is exemplified here by the 
approximate density dependence -Fv,so,j(^/) ~ kj 1 . Such a p~ 1//3 -behavior becomes exact in the 
chiral limit m, = as can be deduced by simple mass dimension counting of the dominant iterated 
l7r-exchange diagrams (the basic argument is that M/f^kf has the correct unit of MeVfm 5 ). 



5 Summary and outlook 

In this work we have calculated the nuclear energy density functional £ [p, r, J] relevant for N=Z 
even-even nuclei in the systematic framework of chiral perturbation theory. Our calculation in- 
cludes the l7r-exchange Fock diagram and the iterated l7r-exchange Hartree and Fock diagrams. 
These few leading order contributions in the small momentum expansion give already a very good 
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equation of state of isospin symmetric infinite nuclear matter ||, |Hj. The step to inhomogeneous 
many-nucleon systems is done with the help of the density-matrix expansion of Negele and Vau- 



therin [ll]]. Our results for the strength functions F T (kf), F^(kf), F so (kf) and Fj(kf) (density 
dependent generalizations of combinations of Skyrme force parameters) are parameterfree. 

We find that the effective nucleon mass M*(p) deviates at most by ±15% from its free space 
value M, with 0.89M < M*(p) < M for p < 0.11 fm~ 3 and M*(p) > M for higher densities 
p < po = 0.174 fm~ 3 . The latter enhancement can be understood from Galitskii's second order 



calculation [13, 14]. Interestingly, a recent large scale fit of (almost two thousand) nuclide masses 
by Pearson et al. \IE\ finds a similarly enhanced effective nucleon mass: M*(po) = 1.05M. 

The strength of the (Vp) 2 -term, F^(kfo), is comparable to that of phenomenological Skyrme 
forces. The magnitude of Fj(kf ) accompanying the squared spin-orbit density J 2 comes out 
larger. 

The strength of the nuclear spin-orbit interaction, F so (kf ), as given by iterated l7r-exchange 
is about half as large as the corresponding empirical value ~ 90MeVfm 5 , however, with the 
wrong negative sign. Since the (positive) contribution from irreducible 27r-exchange to F so (kf) 
is relatively small, there remains the challenge of understanding the microscopic origin of the 
nuclear spin-orbit interaction. Evidently, chiral perturbation theory cannot properly account for 
the underlying mechanisms, whereas relativistic scalar-vector mean-field models give a successful 
phenomenology of the nuclear spin-orbit force. Lorentz scalar and vector mean- fields with their 
in-medium behavior governed by QCD sum rules could also provide the appropriate framework 



for that [IS 



The novel density dependencies of M*(p) and i 7 v,so,j(^/) as predicted by our parameterfree 
calculation should be explored and examined in future nuclear structure calculations (after adding 
a suitable positive constant to F so (kf) in order to minimally repair the spin-orbit coupling). Of 
course one should keep in mind that the prominent low-density behavior of F^{kf) as well as 
Fj(kf) carries little weight in the tails of nuclear density distributions. 

For an extension to even-even nuclei with N > Z the first obvious step would be to include the 
density dependent asymmetry energy A(kf) [§, [K| (subtracted by its kinetic energy contribution) 
in the nuclear energy density functional: 



&as\Ppi Pm 7~p) T n , Jp 1 Jn\ &\Pi ^"i J ] &coul\Ppi Tpi Jp] 
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K } K f OTK f \ 



X ( A( ^ 6M + 12M 3 5QpM 3 j 



(38) 



with p = p p + p n = 2/c 3 /37T 2 , t = Tp + T n and J = J p + J n . In an ordering scheme where 
one counts deviations from homogeneity and deviations from isospin symmetry simultaneously 
as small the energy density functional in eq.(38) would already be complete. However, such a 
formal consideration may be too simplistic in view of neutron skins, neutron halos etc. In any 
case, the density-matrix expansion in eq.(l) can be straightforwardly generalized to the isospin 
asymmetric situation and this way the strength functions of terms like [r n — T p +kj(p p — p n )](p n — P P ), 
(Vp n — Vp p ) 2 , (V Pn — V p p ) ■ (J n — J p ) and (J n — J P ) 2 in the nuclear energy density functional become 
also accessible in our diagrammatic framework. For the Coulomb energy density S cou i[p p ,Tp, J p ] of 
the protons with improved treatment of the exchange (Fock) term, see ref . p0| . 
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